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q< i ABSTRACT 

H ' We extend and reconcile recent work on turbulence and particle heat- 

ing in advection-dominated accretion flows. For approximately equipartition 

magnetic fields, the turbulence primarily heats the electrons. For weaker 

magnetic fields, the protons are primarily heated. The division between 

electron and proton heating occurs between (3 ~ 5 and (3 ~ 100 (J3 = ratio 

>■ ' of gas to magnetic pressure), depending on unknown details of how Alfven 

waves are converted into whistlers on scales of the proton Larmor radius. We 

also discuss the possibility that magnetic reconnection could be a significant 

O ■ source of electron heating. 

OO ' 
On 

O ' 1. Introduction 

Phenomenological models of accretion flows should predict a radiation spectrum for 
a central object of a given mass using a few relevant dimensionless parameters (e.g., the 
^ ■ mass accretion rate in Eddington units, the Shakura-Sunyaev viscosity parameter, a, 

and the ratio of the gas to the magnetic pressure, (3). In advection-dominated accretion 
flows (ADAFs), which are hot collisionless plasmas that have been argued to form around 
some accreting black holes, an additional dimensionless parameter, the fraction of the 
turbulent energy which heats the electrons (= 5), becomes important. In this paper, we 
discuss the electron heating mechanisms in ADAFs; in particular, we show that 5 should 
grow as (3 decreases. 

In ADAFs (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 1995; Abramowicz 
et al. 1995), the accreting gas is unable to cool efficiently and most of the energy 
generated by viscous stresses is stored as thermal energy of the gas and is advected 
onto the central object. As a result, the gas heats up to nearly virial temperatures and 
adopts a two-temperature configuration, with the protons significantly hotter than the 
radiating electrons (Rees et al. 1982; Narayan & Yi 1995). 
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In the context of ADAFs, the issue of which particles (electrons or protons) receive 
the viscous energy acquires particular importance. In "standard" ADAF models (see 
Narayan et al. 1998 for a review) inefficient cooling of the gas occurs because the 
viscously generated energy is assumed to primarily heat the protons. Since the accretion 
flow is effectively collisionless (e.g., Rees et al. 1982) only a small fraction of this energy 
is transferred to the electrons via Coulomb collisions; the total energy radiated by the 
gas (almost all by the electrons) is therefore much less than the total energy generated by 
viscosity. If viscosity were to predominantly heat the electrons, an accretion flow could 
be advection-dominated only at very low densities (and thus very low accretion rates) 
when the electrons themselves are unable to cool efficiently As a result, the optically 
thin ADAF formalism would probably not be relevant for observed systems. 

Recently there has been some theoretical progress in addressing the issue of par- 
ticle heating in ADAFs. The basic physical picture is that ADAFs, like thin disks, 
are magnetized turbulent plasmas. On small scales (much less than the outer scale of 
the turbulence, which is ~ the local radius in the accretion flow), the turbulence will 
be approximately incompressible. Since the Alfven wave is the incompressible MHD 
mode, turbulence on relatively small scales can be described as a spectrum of nonlin- 
early interacting Alfven waves. These nonlinear interactions transfer energy from large 
scale perturbations (waves) to smaller scales until it is dissipated. The relative heat- 
ing of protons and electrons is determined by which particle species is responsible for 
dissipating the turbulence. 

There are (at least) two physical processes which are neglected by simply equating 
particle heating in ADAFs with particle heating by Alfvenic turbulence. The first (which, 
in our view, is potentially important) is magnetic reconnection, which may occur at 
energetically important levels in the accretion flow. This is discussed in §5 and in 
Appendix B. 

The second (which, in our view, is less important) is that, on large scales, the 
turbulence may have a compressible component. This possibility has been considered 
(implicitly) by Blackman (1998), who argues that Fermi acceleration by relatively large 
scale (much larger than p p = the Larmor radius of thermal protons) magnetic pertur- 
bations associated with MHD turbulence preferentially heats the protons to the degree 
required by ADAF models. This conclusion is, however, equivalent to considering the 
collisionless dissipation of the fast mode component of MHD turbulence (which is com- 
pressible), and thus does not apply to the entire turbulent cascade (Achterberg 1981). 
The incompressible component of MHD turbulence, i.e., the Alfven waves, is undamped 
on spatial scales much larger than p p (see §2). 

Two qualitative considerations suggest to us that compressible turbulence may be 
energetically less important than Alfvenic turbulence: (1) the Balbus-Hawley instabil- 
ity, which is likely relevant for the generation of magnetic fields and turbulence in the 
accretion flow, is non-compressive to linear order (Balbus & Hawley 1991). (2) if the 
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turbulent velocity is subsonic, the generation of compressive turbulence may be severely 
suppressed. Throughout this paper, we therefore focus exclusively on Alfvenic turbu- 
lence. Our results can, of course, be trivially recast in terms of the fraction (likely ~ 1) 
of the viscously generated energy residing in Alfvenic turbulence. 

Gruzinov (1998; hereafter G98) and Quataert (1998; hereafter Q98) analyzed par- 
ticle heating by the Alfvenic component of MHD turbulence. In this paper we provide 
a synthesis and extension of their work. We first review Alfvenic turbulence and wave 
dissipation on length scales comparable to or larger than p p (§2). The primary differ- 
ence between Gruzinov and Quataert's analysis is that Quataert's was restricted to these 
"large" length scales, while Gruzinov included a preliminary analysis of what happens 
on smaller scales. In §3 (and Appendix A) we provide a more detailed discussion of 
turbulence and wave dissipation on length scales smaller than p p . We then present a 
model for how the turbulent energy cascades past p p (§4). The principle results of this 
paper are given in §5. This contains our best estimate, along with what we take to be 
plausible uncertainties, of the electron heating rate in ADAFs. 



2. Kolmogorov-Goldreich-Sridhar Turbulence 

To investigate particle heating by turbulence, we must understand both the dy- 
namics of the turbulent cascade, as well as the dissipation mechanisms which operate 
on it. Recent work by Goldreich & Sridhar (1995; hereafter GS) provides the necessary 
characteristics of Alfvenic turbulence. In the MHD limit, linear Alfven waves satisfy the 
dispersion relation u = VA\k\\\, where u is the mode frequency, va is the Alfven speed 
and fc|| is the component of the wavevector along the mean magnetic field. 

GS argue (and recent numerical calculations confirm; Maron 1998) that Alfvenic 
turbulence naturally evolves into a critically balanced state in which the timescale for 
nonlinear effects to transfer energy from a wavevector ~ k to a wavevector ~ 2k (= the 
cascade time, T c ) is comparable to the linear wave period at that scale, T = 2tt/uj; this 
determines how rapid the dissipation must be to halt the cascade. It also implies that 
the cascade is highly anisotropic, with the energy cascading primarily perpendicular to 
the local magnetic field; the parallel and perpendicular sizes of a wave at any scale are 
correlated, with k\\ ~ k 2 / 3 R~ 1 ^ 3 , where R, the outer scale of the turbulence, is ~ the 
local radius in the accretion flow. Since the dissipation of an Alfven wave depends on 
its direction with respect to the background magnetic field (see below), the path of the 
cascade in wavevector space is crucial. 

In collisionless plasmas, such as those in ADAFs, the wave dissipation mechanisms 
of principle importance are wave-particle resonances (molecular viscosity, thermal con- 
ductivity, electrical resistivity, etc. are entirely unimportant). Resonance occurs when 
the frequency of the wave, in the frame moving with the particle along the field line, is 
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an integer multiple of the particle's cyclotron frequency (Q), 

uj — k\\Vn = nQ, (1) 

where v» is the particle's velocity along the magnetic field. For n = 0, resonance occurs 
when the wave's phase speed along the field line, v = ui/k\\, equals dm. A necessary 
(but not sufficient) condition for strong damping is that v be comparable to the thermal 
speed of the particles, so that there are a large number of resonant particles. 

The n = resonance corresponds to two physically distinct wave-particle interac- 
tions. In Landau damping (LD), particle acceleration is due to the wave's longitudinal 
electric field perturbation (i.e., the usual electrostatic force, E z ). The ratio of the elec- 
tron to the proton heating rates (= P) for a wave (with k±p p ^ 1) damped solely by LD 
is Pld ~ ( m eTp I 'vriyTl) 1 ! 2 , which is ^> 1 for the T p T e plasmas of interest to us. In 
transit-time damping (TTD), the magnetic analogue of LD, the interaction is between 
the particle's effective magnetic moment (p = mv\/2B) and the wave's longitudinal 
magnetic field perturbation, B z (Stix 1992). For a wave (with k±p p ^ 1) damped solely 
by TTD, the protons are preferentially heated: Pttd ~ ( r meT e /m p T p ) 1 ^ 2 <C 1.0 This is 
because in plasmas with T p 3> T e , the protons have the larger magnetic moment and so 
couple better to a wave's magnetic field perturbation. 

The Alfven wave has \v\ = va', in plasmas appropriate to ADAFs (T p ^> T e and 
(3 1, where (3 is the ratio of the gas pressure to the magnetic pressure), v a is comparable 
to the electron and proton thermal speeds and so there are a large number of particles 
available to resonate with the wave. In the MHD limit, however, the Alfven wave has 
both E z = and B z = and so is undamped by linear collisionless effects.^ For larger 
wavevectors, the MHD approximations are less applicable and kinetic theory corrections 
to E e and B z become important, leading to finite dissipation. For the perpendicular 
cascade of Alfven waves due to GS, however, when k±p p ~ 1, u> ~ 
flp and son ^ resonances can be satisfied only by particles with |u|| | ~ \tiVL\va/oj ^> va, 
of which there are a negligible number. The cyclotron resonance is thus unimportant in 
dissipating the turbulent energy in the GS cascade. 

In order to accurately assess the properties of Alfven waves with k± ~ p" 1 3> k\\, 
we have solved the full kinetic theory dispersion relation for linear perturbations to a 
warm plasma. These calculations are utilized throughout this paper, but we refer the 
reader to G98 or Q98 for details. Here, and in Figure 1, we summarize the primary 
results. Alfven waves with k± ~ p p 1 ^> k\\ are damped primarily by the protons (via 
TTD). The dissipation rate, 7, is essentially independent of T p /T e , but is an increasing 



1 In this and the previous expression we have taken v the electron and proton thermal speeds; see 
Q98 for a more general expression. 

2 By contrast, the fast mode has B z ^ in the MHD limit and is thus strongly damped by TTD. 
This is why protons are preferentially heated by compressive MHD turbulence in ADAFs. 
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function of (3 (7/0; ~ 0. 16 (3 1 ^ 2 (k^pp) 2 for (3 3> 1 and A;_|_p p < 1). For plasmas with 
(3 1, Alfvenic turbulence is therefore dissipated on length scales <; p p , with most of 
the turbulent energy heating the protons (this is quantified in §5). 

For (3—1, however, the maximal dissipation rate for an Alfven wave in the GS 
cascade (obtained at k±p p = 1) is 7T « 0.1 (see Figure 1). Since the timescale for 
energy to cascade through the inertial range is T c T, this suggests that, for plasmas 
with equipartition magnetic fields (/3 ~ 1), very little of the turbulent energy is dissipated 
on scales comparable to or greater than p p . Consequently, we must investigate the flux 
of turbulent energy past kj_p p ~ 1, as well as turbulence and wave dissipation on scales 
smaller than p p . 

3. Whistler Turbulence 

Alfven waves only exist for kj_p p ^ 1. For kj_p p ^ 1, the same mode is called the 
whistler. In Appendix A we show that all of the turbulent energy which is not damped 
on scales of k±p p ^ 1 is transformed into whistlers at k±p p <: 1 (in particular, we argue 
that there are no other channels through which the energy can travel). 

The following argument regarding particle heating by whistler turbulence can be 
provided. For k±p p ^ 1, but k± ^> k», whistlers have u <C Q p . Thus, in a mode period, 
a particle undergoes many Larmor orbits. Since k±p p ^ 1, the protons (but not the 
electrons) sample a rapidly varying electro- magnetic field in the course of a Larmor orbit. 
As a result, the protons are "frozen out" and become dynamically unimportant (they 
provide a uniform background of positive charge). In particular, they cannot contribute 
to damping the whistler energy. Whistler energy therefore cascades to smaller length 
scales until it is damped by the electrons. 

The more detailed analysis of whistler turbulence and particle heating given in Ap- 
pendix A confirms this picture. The result hinges crucially on our estimate that whistler 
turbulence maintains k\\ <C k± (recall that the GS cascade gives k\\ ~ k^Pp/R) 1 ^ ~ 
10~ 3 k± at k±p p ~ 1). In this case Larmor circle averaging is efficient and the whistler 
energy is dissipated by the electrons (via Landau damping) on scales of kj_p p ~ 30 (see 
Fig. 1). If whistler turbulence were to reach the proton cyclotron frequency before 
k±p p ~ 10 (which we argue is unlikely), the whistler energy would heat the protons via 
the cyclotron resonance (see Appendix A. 3). 

4. Energy Flux through the Damping Barrier 

We therefore assume that most of the energy dissipated on scales ^ p p heats the 
protons, while the energy that gets past k±p p ~ 1 heats the electrons. Here we give a 
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crude model for calculating the damping of the turbulent energy flux; this provides an 
estimate of the fraction of the energy flux that gets through the proton damping barrier 
at k±p p ~ 1, which is needed to calculate the electron heating rate (5). 

We define the full energy spectrum of the turbulence (normalized to the plasma 
density) as E(k±, k\\) and the perpendicular energy spectrum as E(k±): 

< v 2 >= J dk±dk\\E(k±,k\\) = J dk ± E(k L ), (2) 



where v is the plasma velocity. As we are interested in the flux of energy past k±p p = 1, 
we only need E(kj_) for kj_p p ^ 1, which is given by the Kolmogorov-Goldreich-Sridhar 
spectrum 

E(k ± ) = C x ^kl 5 '\ (3) 

where e = cm 2 /s 3 is the energy flux in wavenumber space and C\ is a dimensionless 
constant of order unity (a Kolmogorov constant). There is essentially no dissipation at 
k±p p <C 1, and so e is constant (= e ) in this region. At larger k±, the energy flux e is 
a monotonically decreasing function of k±. The decrease of the energy flux is described 
by 

^ = -2 < 7 (*jl) > E(k ± ), (4) 
where < y(k±) > is the parallel wave number averaged damping rate, 

< lik±) >= eW) ' (5) 

For the GS cascade, the turbulent energy lies inside the cone k\\ ^ u n i/vA, where 
Uni = C 2 e 1 ^ 3 k 2 / 3 is the nonlinear frequency of the turbulence (at scale k±) and C 2 is 
another dimensionless constant of order unity. We model this using a simple expression 
for E(k±, k\\), namely, 

E(k ± , fc„) = E(k ± )S(k {] - u nl /v A ). (6) 

This corresponds to all of the turbulent energy flowing along the line in k space for 
which the linear mode frequency (u = VA\k\\\) is equal to the nonlinear frequency of the 
turbulence (which captures the key physical result of the GS cascade). Using equations 
dD"@) ^ i s straightforward to show that 



e exp 



-2C J d\ogk ± D(k A 



(7) 



where D(k±) = "f(k±, k\\)/u(k±, k\\) (with k\\ = u} n i/vA) is the dimensionless damping 
rate of the mode (the inverse of the quality factor; see Fig. 1) and C = C\C2- 

Recent numerical simulations confirm the GS picture of Alfvenic turbulence (Maron 
1998). They also give values for the dimensionless constants in the above expressions, 
namely C x = 2.5 ± 0.6 and C 2 = 2.2 ± 0.4, so that C w 6. 
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We are interested in the energy flux past k±p p ~ 1. Since D oc k\, the integral 
in equation ([?[) is dominated by the contribution from k±p p ~ 1, which lies outside the 
region where the MHD simulations are valid. Consequently, C ~ 6 must be taken as a 
crude approximation, subject to significant uncertainties. The basic uncertainty is that 
we do not understand the details of how Alfvenic turbulence is converted into whistler 
turbulence on scales of p p . For example, Alfven waves with k±p p ^ 1 can excite whistlers 
with k±p p ^ 1 by three wave interactions (see Appendix A). However, the effective k±p p 
at which this occurs is uncertain. Consequently, so is the upper limit for the integral in 
equation (0). This uncertainty can be absorbed into an uncertainty in C. In addition, 
given that Alfven waves excite whistlers, the precise timescale on which this occurs is 
uncertain (i.e., is the cascade time at k^p p ~ 1 the same as it is in the MHD regime?). 
Small uncertainties in C translate directly (and exponentially) into large uncertainties in 
the energy flux past k^p p ~ 1 and thus into large uncertainties in the predicted electron 
heating rate. 

It should be noted that equation (||) is only a crude approximation for E(k±, k\\)!^\ 
Given the uncertainty in C, however, alternative expressions do not give significantly 
different results for the energy flux past k±p p ~ 1. 

5. Electron Heating Rate 

As before, we define 5 to be the fraction of the turbulent energy which heats the 
electrons; for the model in this paper, it is given by 5 ~ Pttd + e / e o- Pttd is the 
(generally small) contribution to the electron heating from Alfven wave energy dissipated 
at k±p p ^ 1. e/eo is the fraction of the turbulent energy which cascades past the damping 
barrier at k±p p ~ 1 (essentially all of which heats the electrons). This is given by 
equation (|7]) of the previous section. 

Figure 2 shows our best estimates of 5 as a function of (3 for C = 24, 6, and 1.5 
(from left to right). For small (3, the electron heating is dominated by the energy that 
cascades past k±p p ~ 1. This is independent of the proton to electron temperature ratio 
(since the mode dissipation rate at k±p p ^ 1 is; see §2), but is a strong function of (3. 
For larger (3 the dominant contribution to the electron heating is from the dissipation 
of Alfven waves at k±p p ^ 1 (since almost no energy cascades past k±p p ~ 1). Since 
Pttd ~ (rneTe/rripTp) 1 ^ 2 , this is a strong function of T p /T e , but is independent of (3. In 
Figure 2, we have taken T p /T e = 100, appropriate to accretion near the Schwarzschild 
radius of a black hole (where most of the observed radiation originates). Variations in 
T p /T e vertically shift the value of S when it plateaus (i.e., the high (3 values), but do not 
significantly modify the (3 at which the plateau occurs. 



3 For example, as noted above, the GS cascade actually lies inside the cone fen ^ uj n i/vA, not on the 
line fc|| — uj n i/v A . 
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From Figure 2 we infer that turbulence in ADAFs predominantly heats the protons 
(5 ^ 10 -2 , say) only for (3 larger than some critical value, which lies between ~ 5 and 
~ 100 for the values of C taken in Figure 2. These values encompass what we feel to 
be a reasonable estimate of the uncertainty in how the turbulent energy cascades past 
k±p p ~ 1. The corresponding uncertainty in S is extremely large, but this is an accurate 
reflection of the (exponential) sensitivity of our results to unknown details of turbulence 
on scales of the proton Larmor radius. 

The conclusions to be drawn from Figure 2 depend, of course, on the degree of 
one's confidence that all possible sources of electron heating have been accounted for. 
The most accurate interpretation of Figure 2 is that it represents the fraction of the 
small scale (Alfvenic) turbulent energy which heats the electrons. In §1, we argued 
that this interpretation can be broadened to be (within a factor of few) the fraction of 
both compressible and incompressible turbulent energy which heats the electrons (since 
compressible turbulence is unlikely to energetically dominate Alfvenic turbulence and, 
in any case, does not significantly heat the electrons). 

In Appendix B, however, we argue that, at large (3, the physical picture and cal- 
culations leading to Figure 2 neglect a source of electron heating (possibly the most 
important one), namely magnetic reconnection. The reasoning is as follows. Proton 
damping of the turbulent energy at large (3 is essentially a viscous dissipation mecha- 
nism (since the protons carry the momentum of the plasma). For a nontrivial topology, 
the magnetic free energy of the turbulence cannot be fully damped by viscosity. This 
is because changes in the magnetic topology, magnetic flux, magnetic helicity, etc., can 
only be due to resistive effects (regardless of how small the resistivity is). This suggests 
that coincident with the dissipation of turbulent energy by Alfven wave damping is the 
formation of current sheets in which the topology of the magnetic field changes. Crude 
estimates suggest that the energy dissipation in current sheets may be ~ that of Alfven 
wave damping. Unfortunately, we cannot assess whether ~ corresponds to 1 or 10~ 2 , 
etc., and thus we cannot assess whether reconnection is energetically important (that it 
is required for topological changes is clear). One way to address this question is through 
numerical simulations of MHD turbulence with varying magnetic Prandtl numbers (ratio 
of viscosity to resistivity). The relative contributions of Joule heating and viscous heat- 
ing in such simulations could potentially provide important information on the energetic 
importance of magnetic reconnection in large magnetic Prandtl number turbulence. 

6. Conclusions 

Spectral models of ADAFs generally assume that (see, e.g., Narayan et al. 1998) 
(1) magnetic fields are amplified until they are in strict equipartition with gas pressure 
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((3 = 1)[] and that (2) the energy generated by viscous stresses predominantly heats the 
protons (5 & 1CT 2 is a typical value). 

We suggest that these assumptions may be incompatible. Based on an analysis of 
incompressible turbulence (Alfvenic and whistler) and collisionless wave dissipation, we 
find predominantly proton heating only for f3 greater than some critical value, which 
lies between ~ 5 and ~ 100 depending on unknown details of how Alfven waves are 
converted into whistlers on scales of the proton Larmor radius. This does not, of course, 
imply that ADAF models are untenable. Rather, their sensitivity to changes in input 
microphysics (e.g., (3 and 5) should be carefully assessed. This will be pursued in a 
future paper. 

We cannot overemphasize the uncertainty in the numerical values given in this paper 
(which are, e.g., based on the assumption of a uniform thermal plasma). Nonetheless, 
we believe that the basic physical picture of particle heating (wave damping + recon- 
nection), and the general conclusions drawn from it (5 increases as (3 decreases), are 
essentially correct. One point which clearly requires additional investigation is the en- 
ergetic importance of magnetic reconnection (§5 and Appendix B). 

Acknowledgments. EQ thanks George Field, Charles Gammie, Jason Maron, and Ramesh 
Narayan for useful discussions. EQ was supported by an NSF Graduate Research Fel- 
lowship and, in part, by NSF Grant AST 9423209. AG thanks Jason Maron, Peter 
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4 Narayan et al. define the magnetic pressure as B 2 /2Att while the (3 used in this paper defines 
magnetic pressure as B 2 /8n. Thus equipartition for Narayan et al. actually corresponds to (3 = 1/3. 
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A. Whistler Turbulence 

The purpose of this appendix is three- fold. First, we argue that Alfvenic turbulence 
is converted into whistler turbulence on scales of the proton Larmor radius. While the 
details of this process are uncertain, all of the turbulent energy not damped from Alfven 
waves on scales of ~ p p is transformed into whistler energy on smaller scales (there are no 
other channels through which the energy can travel). Second, we discuss the turbulent 
cascade of whistlers. Finally, we discuss the dissipation of whistler turbulence. This, 
together with §2 on Alfvenic turbulence, provides a picture (if approximate and uncertain 
in detail) of turbulence and particle heating from the outer to the inner scale. 

A.l. Alfven — > Whistler: kj_p p ~ 1 

For /3 ~ 1, the full kinetic theory dispersion relation evolves continuously from 
Alfven waves to whistlers as we pass k±p p ~ 1. This is because whistlers are the natural 
generalization of Alfven waves once protons drop out of the small-scale dynamics due to 
Larmor circle averaging. Consequently, Alfven waves at k±p p ^ 1 can excite whistlers 
at k±p p ^ 1 by three- wave interactions (i.e., the resonance conditions can be satisfied) .0 
For (3 ^ 100, the situation is more complicated because there is a region (the "damping 
barrier") of wavevector space (0.85 ^ k±p p ^ 1.15 for (3 ~ 100 and larger for larger 0) 
in which Alfven waves do not propagate (see G98). Because this region is rather narrow 
for the 13 of interest to us, Alfven waves can excite whistlers across the damping barrier 
(narrow meaning that the jump in k± is less than a factor of ~ 2). This establishes that 
the turbulent energy not damped on scales of k±p p ~ 1 can be converted to whistler 
energy on smaller scales (although, as discussed in §4, the precise details - e.g., the 
timescale - are uncertain). 

To establish that no other channels of energy travel are possible is straightforward. 
Whistlers are the only modes with k±p p ^ 1 and uj <C Q p ("sound waves" are too 
strongly damped to be excited in collisionless plasmas; see Q98). They are therefore the 
only possible sink of the Alfvenic energy not damped on scales of k±p p ~ 1. 

A. 2. Whistler Cascade 

Like collisionless Alfven waves, collisionless whistlers can be described hydrody- 
namically. Hydrodynamic equations valid at kj_p p ^ 1 can be obtained as follows. The 
protons are dynamically frozen; their sole function is to create a positive charge back- 



5 Note that because the GS cascade is strong, the frequency "resonance" condition is actually quite 
broad. 
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ground on which the electrons and magnetic fields evolve. Electrons move freely along 
magnetic field lines and so En = 0. In the perpendicular direction, electrons move with 
the ExB-drift velocity, vi = cE x B/B 2 . These two equations give 



E + -vxB = 0. (Al) 

c 



Neglecting displacement currents, 



47T 

V x B = nev. (A2) 

c 



From equations (|A1|) and flS~2), we obtain the Electron Magnetohydrodynamics (EMHD) 



equation (Kingsep, Chukbar, & Yan'kov 1990) 



d t B = V x (B x V x B). (A3) 

Anne 



The linear waves in equation (S3) are whistlers; their dispersion relation can be presented 

as 



v 



2 



\(k±p P )k\\\, (A4) 

where v p is the proton thermal speed. 



v p 



To see at what scales the protons freeze out, we compared the whistler dispersion 
relation obtained numerically from a full plasma permittivity tensor (including both 
electrons and protons) to the analytical dispersion relation flS3|). At (3 = 100, the 
dispersion relations agree to 30 % at k±p p = 1.5, to 10 % at kj_p p = 2, and to 3 % at 
k±p P = 3.Q 



The turbulent cascade in (A3) was described by Kingsep, Chukbar, & Yan'kov 
(1990). A very clear discussion is given by Goldreich & Reisenegger (1992). These 
authors assumed that whistler turbulence is roughly isotropic in wavenumber space, and 
that the turbulence is weak. Their argument for the turbulence being weak is analogous 
to the standard argument for the "weakness" of Alfvenic turbulence. Assuming strong 
isotropic turbulence, one calculates a nonlinear frequency at a given scale. This frequency 
turns out to be smaller than the linear wave frequency; the turbulence must therefore 
be weak. 

As explained by GS for the MHD case, however, the above "explanation" does 
not necessarily work. Instead of being weak and isotropic, the turbulence is strong 
in a narrow cone in wavenumber space (k» <C k±), for which the linear frequency is 
smaller than the nonlinear frequency. Numerical simulations seem to confirm the GS 
picture in the Alfven wave case (Maron 1998). But the same picture must hold true 



6 Other relevant plasma parameters in this example are v e /c = 0.5, v p /c — 0.33, fen <C k±. For details 
see Q98 or G98. 
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(to some extent) in the whistler case. In a narrow cone in wavenumber space, the 
turbulent cascade can proceed at a (fast) nonlinear rate. Outside of the cone, the cascade 
slows down (the turbulence is weak). This is because nonlinear interactions outside the 
cone must satisfy the frequency resonance condition (interactions are restricted to the 
resonant manifold). 

For our problem, energy is injected into whistler turbulence in the narrow cone for 
which the turbulence is strong (because the Alfvenic turbulence which excites whistler 
turbulence is itself strong). This suggests that a large part of the energy will continue 
to cascade in the strong cone, maintaining k\\ <C k±. Numerical simulations in two di- 
mensions (Biskamp, Schwarz & Drake 1996) confirm that strong Kolmogorov turbulence 
results in both the Alfven and the whistler cases. 

We have argued that the energy in whistler turbulence injected at k\\ <C A;_l preserves 
this anisotropy because the turbulent energy cascade is faster in the strong cone where 
the nonlinear frequency shifts exceed the linear frequency. We cannot, however, rule 
out that the strong cone may be leaky. It is a peculiar feature of Alfvenic turbulence, 
not shared by its whistler counterpart, that the 3-wave (and 4- wave) resonance condi- 
tions drive the turbulent energy into the strong cone, regardless of their starting point 
in wavevector space. Our analysis of whistler turbulence is therefore only suggestive. 
Thankfully, the precise path of whistler turbulence in wavevector space is not needed to 
assess its particle heating properties. All that is needed is that it maintains k\\ <C k± 
(this is quantified below). 



In this subsection, we assess the dissipation of whistler turbulence. Assuming that 
whistler turbulence maintains kn <C k±, we show that the protons cannot be heated by 
whistlers. 

For a (subthermal) wave with k±p p ^ 1 damped solely by transit time damping, 
the relative heating of electrons and protons is given by 



The (kj_p p ) 3 factor describes quantitatively the effects of Larmor circle averaging (com- 
pare with the k±p p & 1 expression in §2). For T p ~ 100T e , electrons are heated more 
than protons for k±p p ^ 7.f\ The decrease in the whistler dissipation rate for k±p p <; 1 
in Figure 1 is due to the freezing out of the protons. Even though the protons are in 
principle heated more than the electrons for k^p p ^ 7, the decrease in the dissipation 



7 Note that for all k±p p Landau damping preferentially heats the electrons for T p ^> T e . 
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rate ensures that very little energy is dissipated between 1 & k±p p ^ 7. For length scales 
smaller than kj_p p ~ 7, Larmor circle averaging entails that protons cannot be efficiently 
heated by the n = Landau resonance (see eq. |]T|). This conclusion is essentially 
independent of the details of the whistler cascade. 

For k^p p ^ 1, protons can also be heated by the cyclotron resonance (for wave 
frequencies ~ the proton cyclotron frequency; the \n\ = 1 resonance in eq. |l|). For 
large k±p p , the dominant electron heating mechanism is Landau damping. The relative 
importance of proton heating by the cyclotron resonance and electron heating by Landau 
damping is given roughly by (m p /m e )(k±p p )~ 3 , where we have taken v a ~ v p ~ v e ~ 
v ~ c and to ~ 

For k±p p ^ 10, Larmor circle averaging is sufficiently strong that proton heating by 
the cyclotron resonance is less important than electron heating by Landau damping.^ 
Only if whistler turbulence reaches uj ~ Q p before k±p p ~ 10 can the protons be heated 
by whistler turbulence. This would require a complete reversal in the direction of the 
turbulent cascade since at k±p p ~ 1 the GS cascade yields uj ~ 10~ 3 fi p . Our estimates 
of whistler turbulence in the previous subsection suggest that it reaches uj ~ Q p only 
for k±p p ^> 10. Consequently, Larmor circle averaging of the protons is efficient and the 
whistler energy is damped by the electrons (by Landau damping) on scales of k±p p ~ 30 
(see Figure 1). 



B. Reconnection 

In this appendix, we discuss why reconnection (Joule heating) may be the primary 
heating mechanism for the electrons, especially at low magnetic fields (when whistlers are 
not significantly excited). We argue that, to order of magnitude, comparable amounts 
of energy are damped by viscous heating of the protons and by Joule heating of the 
electrons. Bisnovatyi-Kogan & Lovelace (1997) first suggested that Joule heating of the 
electrons could be important in ADAFs (we do not, however, agree with the particular 
scenario they propose; see Q98 or Blackman (1998)). Blandford (1998) also suggested 
that reconnection could be important. The following analysis attempts to assess the 
potential significance of reconnection in a semi-rigorous manner. 

For (3 3> 1, Alfven waves are damped at k±p p < 1 and the plasma dynamics can be 
described by the MHD equations. There are two heating mechanisms in the framework 
of MHD - viscous heating and Joule heating. Roughly speaking, viscosity (u) heats 
the particles that carry the momentum, i.e., the protons, and resistivity (rj) heats the 



8 In deriving this expression we have used an analytical approximation for the relative strengths of 
the perpendicular and parallel electric field components of the whistler wave (from G98). 

9 Our numerical calculations with the full plasma permittivity tensor confirm this result. 
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particles that carry the current, i.e., the electrons. We know that j/ > ?] in a /3 > 1 
plasma - protons are heated when Alfven waves are damped by viscosity. 

Is it true, however, that Joule (electron) heating dissipates much less energy than 
viscous (proton) heating in a v rj plasma (i.e., for large magnetic Prandtl numbers)? 
In other words, can Alfven wave damping actually dissipate all of the turbulent energy? 

The magnetic field of the accreting plasma is ostensibly self-generated. The mag- 
netic field energy is constantly being created, destroyed, and created again. To demon- 
strate that reconnection is unavoidable, consider the relaxation of a magnetic field per- 
turbation created "by hand" at time zero. This initial magnetic field is assumed to be 
a generic solenoidal field with a characteristic strength B and a characteristic scale L; 
we also assume that the normal component of the field is specified at the boundary of a 
box with side ~ L. Now let the field evolve according to MHD equations with v ^> r\. 

Magnetic pressure and tension cause plasma motions on scales ~ L. A Kolmogorov- 
GS cascade develops and magnetic energy is converted into kinetic energy; viscous heat- 
ing by the protons dissipates the kinetic energy. In a time ~ L/va, the kinetic energy 
goes to zero and the magnetic field reaches an equilibrium state defined by 

V x B x B = AnVp, (Bl) 

where p is the plasma pressure. Generically B 2 ~ B 2 ~ (B — B) 2 so that the energy 
heating the protons is ~ B 2 L 3 . 

As explained by Arnold (1986), the equilibrium magnetic field lines described by 
equation (Bl) have a very special topology. They lie on a set of nested tori with p = 
const. For rj — 0, however, the magnetic field is frozen into the plasma (magnetic flux, 
helicity, and the topology of the field lines are all conserved). Since the initial magnetic 
field topology was a generic one (by assumption), it is impossible to reach the equilibrium 
configuration with rj = (since that requires changing the topology of the field). On the 
other hand, the equilibrium must be reached, because out of equilibrium the magnetic 
field causes plasma motions and therefore viscous heating. 

The resolution of the paradox is the formation of current sheets. Topological barriers 
on the way to equilibrium are squeezed to zero volume. No matter how small, resistivity 
becomes important in the current sheets thus allowing reconnection (= changes in the 
magnetic field topology). Generically the magnetic free energy released at this stage is 
~ B 2 L 3 . 

The only way to avoid reconnection is to assume that the magnetic field has a 
trivial topology. For example, one can assume that magnetic field lines constantly lie 
on surfaces that are topologically equivalent to a set of nested tori. The dynamo nature 
of the magnetic field in accretion flows, however, seems incompatible with a trivial 
topology (cf the complex structures observed in numerical simulations of the Balbus- 
Hawley instability). 
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While we do not know how to rigorously estimate the fraction of the turbulent 
energy damped in reconnection events, the potential importance of reconnection for 
electron heating, and its unavoidability even in large magnetic Prandtl number plasmas, 
seems clear. 
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Fig. 1. — Dimensionless dissipation rate of the Alfven (k±p p < 1) and whistler (k±p p > 1) 
modes, taking k\\ k±. The solid (dashed) curve is for /3 = 1(10), while both curves 
take T p = 100T e . The peak in the dissipation at k±p p ~ 1 corresponds to proton heating 
while the strong damping at k±p p > 30 corresponds to electron heating. 
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Fig. 2. — Estimates of 5, the fraction of the turbulent energy which heats the electrons, 
versus f3 (for T p = 100T e ). The three curves correspond to (from left to right) C = 24, 
6, and 1.5 (see eq. M). 



